*Creates Figures 1a and 1b from BMSVY(2020)
version 16.0
u figure1,clear
graph twoway (scatter change_ltfp_w change_management if treatment==1,ms(x)) (scatter change_ltfp_w change_management if treatment==0,ms(o) ),  xlab(-0.1(0.1)0.7) ylab(-1(0.5)1.5) graphregion(color(white)) legend(off) saving(Figure1a,replace)
graph twoway (scatter change_lqdi_w change_management if treatment==1,ms(x)) (scatter change_lqdi_w change_management if treatment==0,ms(o) ),  xlab(-0.1(0.1)0.7) ylab(-2(0.5)1) graphregion(color(white)) legend(off) saving(Figure1b,replace)
	